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■ Abstract 

We study Kraichnan's model of a turbulent scalar, passively advected by a Gaussian 

. random velocity field delta-correlated in time, for every space dimension d > 2 and eddy- 

O I diffusivity (Richardson) exponent < C < 2. We prove that at zero molecular diffusivity, 

■ or K = 0, there exist unique weak solutions in (£7®^) to the singular-elliptic, linear 



PDE's for the stationary A^-point statistical correlation functions, when the scalar field 



' is confined to a bounded domain with Dirichlet b.c. Under those conditions we prove 

5t ! that the A^-body elliptic operators in the spaces have purely discrete, positive spectrum 

and a minimum eigenvalue of order L^'^, with 7 = 2 — <^ and with L the diameter of fl. 
We also prove that the weak L^-limits of the stationary solutions for positive, pth-order 
hyperdiffusivities Hp > 0, p > 1, exist when Hp and coincide with the unique zero- 
diffusivity solutions. These results follow from a lower estimate on the minimum eigenvalue 
of the iV-particle eddy-diffusivity matrix, which is conjectured for general N and proved in 
detail for = 2, 3, 4. Some additional issues are discussed: (1) Holder regularity of the 
solutions; (2) the reconstruction of an invariant probability measure on scalar fields from the 
set of A^-point correlation functions, and (3) time-dependent weak solutions to the PDE's 
for A^-point correlation functions with initial data. 
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1 Introduction 

We study the model problem of a scalar field 9{r, t) satisfying an advection-diffusion equation 

{dt + y^-Vr)e = KAr9 + f (1) 

in a bounded domain Q of Euclidean d-dimensional space R'^, with Dirichlet conditions on the 
boundary d^l. The scalar source /(r, t) is assumed a Gaussian random field, white-noise in time 
but regular in space. Precisely, we take / with mean (/(r,t)) = /(r) G L^(r2) and covariance 

{fir,t)f{r',t')) - {f{r,t)){f{r',t'))=Fir,r')6{t-t') (2) 

with F £ (J7 $7). The velocity field is also assumed Gaussian, white-noise in time, zero- 
mean with covariance 

{v^{r,t)vj{r',t')) = Vij{r - r')6{t - t') (3) 

The velocity to be considered is a divergence-free random field in R'^ and, for convenience, 
statistically homogeneous. There is no reason to insist on Dirichlet b.c. for the velocity field. 
The spatial covariance matrix V we consider is defined by the Fourier integral 



V.i(r) = Do -—, Ik'+m^) PMk)e-''-. (4) 



where < C < 2 and P^ili.) is the projection in R"' onto the subspace perpendicular to 
k. This automatically defines a suitable positive-definite, symmetric matrix-valued function, 
divergence-free in each index. The model originates in the 1968 work of R. H. Kraichnan |l|] 
and has been the subject of recent analytical investigations [§, |3|, ^ ^, ^, 0, It is not hard to 
show that 

+ ■■■ (5) 



asymptotically for mr <^ 1, with Vq and Di constants proportional to Dq, given below. See 
also Section 4.1 of Q. The exponent (" has the physical interpretation of an "eddy-diffusivity 
exponent" analogous to the Richardson exponent 4/3 M. 



The remarkable feature of Kraichnan's model, which makes it, in a certain sense, "exactly 
soluble" is that A^-th order correlation functions 0Ar(ri, r^r; t) = {6{ri,t) ■ ■ ■ 9{r]\f,t)) satisfy 
closed equations of the form 

dtQN = -h'n^&n + ^7(r„)ejv-i(...f;...) 

n 

n ) I'm ■ ■ ■ I'm • • • )• (6) 

pairs {nm} 

In this equation for the A^-correlator only itself and lower-order correlators appear ^, |3|, |^. 
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Here Tl/^ is an elliptic partial-differential operator in fi®^ defined as 

d N 



Vij{rn 



'■mj 



N 

^E^r„, (7) 



n=l 



with Dirichlet b.c, where Xin are Cartesian coordinates in (R ) • However, the operator 
obtained by taking k — > is degenerate, i.e. it is singular-elliptic. We refer to "Hat as the 
A^-body convective operator because it accounts for the effects of the velocity advection alone 
in the equation (^) for A^-point correlations. Because of the degeneracy for k ^ 0, the solutions 
of the parabolic equation are expected in that limit to lie only in a Holder class with 
7 = 2 — C. As the differential operator is of second-order, these solutions must then be taken in 
a suitable weak sense. Despite the degeneracy, the linear operator 7i]\f is formally self-adjoint 
and nonnegative in the inner product of functions on 0®^. This suggests that an L^-theory 
of weak solutions to Eq.(^ may be appropriate. We shall develop here such a theory in detail. 
The key to the analysis of the k ^ limiting solutions is a proof of existence and uniqueness 
directly for n = 0. 

Let us state precisely the main theorems of this work. We shall actually consider a somewhat 
more general model than Eq.(|l]), namely, 

{^t+^r■Vr)e = -Kp{-Arfe + f (8) 

with p > 1, in which Hp is a so-called hyperdiffusivity of order p. This allows us to establish a 
universality result concerning the independence of limits on p. In this case, the closed correlation 



i,j=l ?i,m=l 



AT 



equations (^) are still satisfied, with the operator (0) replaced by 

^^,(rn-r™)^- +Kp^(-Ar„r (9) 

Note that this operator requires higher-order Dirichlet b.c, namely, elements in its domain 
must have zero trace on the boundary for the first k = \p — (1/2)] derivatives. However, our 
first main result is for the solution of that equation directly at Kp = 0: 

Theorem 1 Assume that d > 2 and < C < 2. Then, for integers N > 1, the equation ^ at 
K = has a unique stationary weak solution in L Away from the codimension-d 

set where pairs of points in 11= (ri, ...,rjv) coincide, the solution 0J^(R) is in Hq ( il®^ 



This ideal zero-diffusivity solution is, in fact, the physically relevant one in the limits Kp — > 0, 
as shown by our second main result: 

Theorem 2 Assume that d> 2 and < C < 2, and also p > 1. 

(i) For integers N >1, the equation with TL^'^'> generalized to TL^'^^'\ has a unique stationary 
weak solution Q^^^^* in {o!^^^ , which, in fact, belongs to the Sobolev space Hq ^$7®^^ . 
(a) The weak-L?' limit exists as Kp ^ and w — lim^p^o ©jy''^* ~ ®jV' 

To prove these results requires a spectral analysis of the A^-body convective operator TCn- In 
fact, we show that this operator has pure point spectrum, using a criterion borrowed from a work 



of R. T. Lewis [10|. Discreteness of the spectrum was already shown by Majda |Q] in his simple 
version of the model. For our theorems above, we do not really require that TC^ have a compact 
inverse, but merely a bounded inverse. To prove this, we require an estimate from below on 
the quadratic form associated to TCn- This is proved in two steps. First, for each integer 
> 1 we define the (Nd) x (A^d)-dimensional matrix [GN(R}]in,jm = {vi{Yn)vj{vm)) , i,j = 
1, ...,d,n,m = 1,...,N. Physically, this is interpreted as an N-particle eddy-diffusivity matrix. 
Mathematically, it is the nonnegative Gramian matrix of the Nd elements fi(r„) in the 
inner-product space of the random velocity field. It is nonsingular if and only if these Nd 
elements are linearly independent. We shall prove below (Proposition 2) that its minimum 



eigenvalue obeys A^™(R) > CAr[p(R)]^, where p{K) — min^^m Ir^ — r^l, when N — 2,3,4. 
The second step of the proof uses only this property of GAr(R), which is conjectured to hold 
for all > 1. As a consequence of this estimate, we prove a lower bound on the operator 
quadratic form, reminiscent of the well-known Hardy inequality |11] (Theorem 330). For the 



operator with Dirichlet b.c. we may adapt a convenient proof of the Hardy-type inequality due 
also to Lewis Unfortunately, as explained below, this proof does not work with periodic 
b.c. although the inequality is likely to hold there as well (for zero-mean functions). Lewis' 
argument is also too restrictive to permit treatment of other models with more natural b.c. on 
the velocity field. In a real turbulent flow with velocity field governed by the Navier-Stokes 
equation, the realizations of the velocity field would satisfy also Dirichlet b.c. This behavior 
may be mimicked with the Gaussian random velocity fields by taking as their covariance 

VPir,r') = Af,(r)K,(r - r')An(r'), (10) 

in which An(r) is a suitable "wall-damping function". It should be taken as some decreasing 
function of the distance to the boundary dQ, vanishing there as some power. Of course, with 
this choice of velocity covariance, a lower bound directly follows from our present work that 
A^'°(R) > CAr[p(R)]^[Af7(R)]2, where Aq(R) = mini<„<Ar An(r„). While we expect the main 
results of this work to carry over to such models, it requires a different proof of the generalized 
Hardy inequality. We will return to this problem in a later work. 

Let us summarize the contents of this paper: In Section 2 we establish the required prop- 
erties of the model velocity covariance and the resulting A^-particle eddy-diffusivity matrix, in 
particular the lower bound on the minimum eigenvalue. In Section 3 we study the operator 
quadratic form, and prove its principal properties, such as the generalized Hardy inequality. 
Finally, in Section 4 we exploit these results to prove the main Theorems 1 and 2 above. In 
the conclusion Section 5 we briefly discuss three other problems: regularity of solutions, the 
reconstruction of an invariant measure from the stationary A^-point correlation functions, and 
time-dependent solutions to the parabolic PDE's for the A^-correlators. 



2 Properties of the N-Particle Eddy-DifFusivity Matrix 

(2.1) The Velocity Covariance Matrix 

We first state and prove the regularity properties of the velocity covariance matrix elements 
(Vij(r)) that we will need for later analysis. We have made the choice of Eq.(|^ just for 
specificity. In fact, any velocity covariance with the following properties would suffice. 

Lemma 1 The elements of velocity covariance matrix Vij{r), r G R'^, are C°° in r ifr ^ 0, and 
near r = 0, with £ (0,2). Moreover, there is a positive number po such that if r £ [0, Po]; 
we have the local expansion: 

Vj{r) = Vo6ij-D^-r< ■ S^j + (^dij - + O (m V) (11) 

Proof: The matrix Vij{r) can be written as 

Vij{r) = V{r)5^,+d^djW{r), (12) 

where the function V{r) is defined by the integral 

and W{r) is given by the (for d = 2, principal part) integral 

so that — AW = V . The scalar function V{r) is essentially just the standard Bessel potential 
kernel |jl2|, and may thus be expressed in terms of a modified Bessel function: 

2i-(C/2)™-C 

The Hessian didjW{r) of the function W of magnitude r = |r| alone is 

dA 

didjW{r) = 6ijA{r) + fifj • r—{r), (16) 
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with A{r) = W'{r)/r and r = r/r. However, because Tr (V VM^) = —V, a Cauchy-Euler 
equation follows for A{r): 

fJA 

r—{r)+d-A{r) = -V{r). (17) 
dr 

Due to the rapid decay of its Fourier transform, the function A(r) is continuous. Thus, the 
relevant solution is found to be 

A{r) = -r"^ r p'^~^V{p)dp. (18) 
Jq 

in terms of V{r). Using this expression for A{r), along with Eg. (p!6|) , we thus find 

y,,(r) = {V{r) + A{r))5^j - {V{r) + d ■ A{r))?,rj, (19) 

for Vij as a linear functional of V. If V has a power-law form, V{r) = Br^, then it is easy to 
calculate that 



, ... (20) 

By means of the known Frobenius series expansions for the modified Bessel functions (e.g. |13| 
(9.6.2), (9.6.10)), it follows that 



-^i^.(-) = S-^:^-" + o (21) 



From these terms for K^{z) we obtain, upon substituting Eq.(|l5|) into Eq.(p^), the claimed 
asymptotic expression for Vij{r) in Eq.(|Tl|), with 

Vo = Do ^f-v • m-f, (22) 

and 

(d-i)r(Y) 

Di = Do /, . X • (23) 

(47r)'^/2.2C.C-r(^) 

Finally, the Bessel function Ku{z) is analytic in the complex plane with a branch cut along the 
negative real axis. Thus, the stated smoothness properties of Vij follow. □ 

We shall denote the second term on the right hand side of ( pT] ) as —r^Qij. Obviously, (Qij) is 
positive definite uniformly in r. We will denote by r„m = r„— the vector, and r„m, = |r„— r^l 



the scalar distance from r„ to r^; Vij the matrix elements, and Vnm the matrix evaluated at 
^nm- We show two more lemmas. 

Lemma 2 Let r^, r = 1,2,3, be any three points in R!^, and ri2 < ris, r\2 < ^23- Then there 
is a constant C depending on po and ^ in Lemma 1 but independent of r 12, ris, and r23 such 
that: 

|Vi3 - V23I < C'ri2max(r^^\r|^^). 
Proof. If C e (1, 2), then VV G C^~^, and so by Lemma 1: 

|Vi3 - V23I = |ri2 • VriVlrJ = |ri2 • (VriV|r, - VriV|r=o)|, 

<crv2rl <cri2max(rj2 ,r^3 ), 

where vq = Qv\ + (1 — ^)r2, for some d G (0, 1). The case = 1 is obviously true by the mean 
value theorem. Now if C ^ (Oj 1)) niax('"i3) ''23) > PO) then using V G away from zero, we 
have: 

|Vi3 - V23I < cri2 < cri2(mmax(ri3,r23))''~"^ 
< c{po,m)ri2 max(r^3"\r^3"^). 
If C £ (0, 1), and max(ri3, r23) < po, we employ local expansion to calculate for any x 7^ y: 



+c\y<& - Ml 

< cmax(x^~^, y^~^)|x — y| + cy^ 
The latter term is just: 



2 2 
XiXjy yiyjx 



- c\ i^i^j ~ yiyj)y + viVjiy - ^ ) i 
- nn,c ( \^-y\ , \^-y\y , \^-y\{x + y) 

V kc x-^ 
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With no loss of generality, we assume that y < x; otherwise, we simply switch x and y. It 
follows that 

\Vij{-x.) - Vij{y)\ < c|x - y| max(a;^~\/~^) +c|x-y|a;^~^ < c|x - y| max(x^~\ 
We complete the proof with x = ri3, and y = r23. 

Lemma 3 Assume that r\2 < r^; ri3 = 0(ri4) = 0(r23) = 0(r24); ^ < e G (0,1). Then 
there exist cq and a positive constant ci depending on pg, (, maximum and minimum ratios of 
ri3, ri4, r23; and r24, such that: 



for all e G (0, eo). 

Proof: Applying the mean value theorem to i^(ri) = V13 — V14, we get for = ^ri + (1 — 9)r2 
that: 



By the smoothness of Vn Vij when the distance of ri from rj, z = 3, 4 is larger than ^ (which 
is possible if e is small enough), we obtain: 



Vi3 - Vi4 - (V23 - V24)| < cirururis , 



F(ri) - F(r2) = ■ V,,F 



If max(ri3,r24) > ^, then 



Vri-Flr, = Vi3 - V14 



ri=re- 



< cip^rzA 



from which it follows that: 



F(ri) - F(r2)| < Ciri2r34 < Ciri2r34r|3 



On the other hand, if max(ri3,r24) < we use local expansion in Lemma 1 to get for each 



matrix element: 



(F(ri)-F(r3))i,- = {-D,)riSij + 



d-l 



c 



- 




)] 
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- (1-2) 

= (-^i)ri2-Vr,(4[<5., + ^(5,,-^^; 
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- rl[6,, + -^{6.,-^^mn=re), (24) 

where the notation (1 — > 2) means the same terms as before except that subscript 1 is replaced 
by 2. Let us calculate the ri gradient in (p3) as (fc meaning the fcth component of this gradient): 



Jk) (k) ^ (i) (j) ^ (i) (j) 



/Jk) (k) ^ U) (j) \ 

/- / o («) (i) C^) X U) \ 
+ ^ 4^^™^+4^ + 4i^-(3^4) . (25) 

" \ '13 '13 '13 / 

Note that the first term of the right hand side of ( p5|) is bounded by: 

C{C,d)\4-'-4-'\ < C{C,d)4frs,. 

We can think of 

Jk) _ Jk) ^ (i) (j) 

^ /<: ^13^13 M 



ri3 a - I rf3 

as a bounded function of the unit vector ri3 along ri3. Hence the second term of (p5|) being 



the difference of two values of this function at two points ri3 and ri4 is of the order O(f^) 
Thus the second term is bounded by 

- C-i -1 ^ - C-2 



Similarly, the third term is bounded as such. Combining the above with (24) we deduce that 
l-F(ri) — -F(r3)| < ciri2r'iir\^ . The proof of the lemma is complete. 

10 



(2.2) The N-Point Eddy-Diffusivity (Gramian) Matrix 

As in the Introduction, we define for each integer > 1 the (Nd) x (A(i)-dimensional Gramian 
matrix [GN(R}]in,jm = {vi{'rn)vj{rm)) ■ For the moment we consider general velocity covari- 
ances, given by a Fourier integral 

^..(r) = y V-,(k)e^i'-, (26) 

with V(k) > for each k G R'^. The basic properties are contained in: 

Proposition 1 For each N >2 the matrix Giv(R) has the following properties: 

(i) G^(R) > 0. 

(ii) Assume that for all k G R'^ the velocity spectral matrix V(k) > on the subspace 
orthogonal to the vector k. In that case, Giv(R) has a nontrivial null space if and only if 
Tn = Tm for somc pair of points n ^ m. 

(Hi) For the same hypothesis as (ii), i/ {ri, tat} has K subsets of coinciding points, with 
Nk points in the kth subset, k = 1,...,K, then the dimension of the null space of Gn(R) is 
Y^k=i{^k — ^)d. The null space consists precisely of vectors S = {^i, ...,^^) with the property 
that 

E ^n, = 0, (27) 

nfe=l 

for each k = 1, where the sum runs over the N^. coinciding points in the kth subset. 

Proof, (i) Obvious from the stochastic representation. (ii)Sz (Hi) Let us assume that the Ad- 
dimensional vector H = (^1, ^tv) belongs to KerGAr(R). Then, using the definition of GAr(R) 
and the Fourier integral representation Eg. (|26|) , it follows that 

= (H,G^(R)H) = / ^ (^f: ^„e^k.r„j.v(k). (^fj ^.e'*^-"^ . (28) 

This can only occur if the nonnegative integrand vanishes for a.e. k G R'^. Because of our 
assumption on V(k), this implies that 

N 

5:^„e^'^-- = a(k).k, (29) 

n=l 
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for a.e. k G R'^ with some complex coefficient a(k). Taking the vector cross product with 
respect to k and then Fourier transforming, we obtain that 

N 

^^„xV<5(r-rJ = 0, (30) 

n=l 

in the sense of distributions. Therefore, for any smooth test function 99, 

TV 

5]^„x(V¥.)(rJ = 0. (31) 

n=l 

Because the values of V(/3 may be arbitrarily specified at any set of distinct points, it follows 
that 

K ( Nk \ 

E E^nJxa'^ = (32) 

k=l \nfc=l / 

with afc G R'' arbitrary. This immediately implies that Eq.(27) is both necessary and sufficient 
for H to belong to KerG7v(R). Furthermore, this subspace has dimension '^k=ii^k — ^)di which 
completes the proof of (in). Finally, (ii) follows from (in) by observing that KerG7v(R) = if 
and only if K = iV and iVfc = 1 for all k = I, N. □ 

For the particular choice of covariance function defined by Eq.(^) for < C < 2, we need 
also the following crucial lower bound: 

Proposition 2 For each < C < 2 and d > 2, there exists for each N > 2 a constant 
Cn = C]\f{d,(^) > so that the minimum eigenvalue Ajy™(R) o/G7v(R) satisfies 

Ar (R) > Cn ■ [p(R)]^, (33) 

with p{K) = mXD.n^rnrnm- 

The above property will be proved in detail in this paper for = 2, 3, 4. While the proof in 
these cases strongly suggests the result is true for all > 2, the argument becomes increas- 
ingly complicated for larger values of A^. We shall leave the discussion of the general to a 
future publication, although we point out that many parts of the argument below apply for 
the general case. Note that we can view Gat as a matrix parametrized by the C, power of the 
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minimum distance, e = p^. Let A™™ = \n{^) be the minimum positive eigenvalue of Gat with 
corresponding unit eigenvector H™™ = Hiv(e). Then, by the standard formulae of degenerate 
first-order perturbation theory (see Kato, p^): 

Ajv(e) = (H^(0), G^(e)3^(0)) + ©(e^). (34) 

We have used the fact that AAr(e) is at least twice differentiable in e near zero: see Theorems 
II. 1.8 and II. 6. 8. Furthermore, H7v(0) is in the null space of GAr(O). Thus, by Proposition 
l(iii), H7v(0) = (^i,"-,^Ar) such that J2n=i^n — 0- By simply minimizing over this entire 
subspace of vectors H, we shall show that the righthand side quadratic form of (|3^, denoted 
by Qn{$,i-, ■ ■ ■ ,^Ar_i), is bounded from below by a constant times e. Thus A(e) obeys the same 
type of lower bound. 
Proposition 2, N=3 Case 

Remark: The following proof for Proposition 2, = 3, also implies the lower bound C2f'i2 
for the N = 2 case. 

Proof: Let r„, n = 1,2,3, be three distinct points in R"', d>2. Then we show that there is a 
positive constant C3 = 6*3(^0)) where po is the scale of local approximation (p!l[), such that the 
minimum eigenvalue of G3 is bounded from below by C^p^ . It suffices to treat the situation 
where p < po, otherwise, we conclude with Proposition 1. Let Cq be a large but 0(1) constant 
to be determined, and let ri2 = p for definiteness. 

Case I: Suppose now that ^ < Cq, and ^ < Cq. By further reducing the size of p, we can 



ensure that pCo < pQ. Now write: 

V -^1-^2 J 

then: 



1 
1 
-1 -1 1 



^2 

V oy 



Q3 = ((^i,e 



2(V(0)-Vi3) V(0)+Vi2-Vi3-V 



2^ 



V(0) + V12 - Vi3 - V: 



23 



2(V(0) - V23 




13 



Since all the three distances are less than po, we apply lemma 1 to see that ^^^^^^ ^'■^^ < Cq. 
Therefore we can factor out . The remaining entries are bounded by Co, and we also know 
that they form a positive definite matrix. Hence by continuity of eigenvalues on the matrix 
entries, we get the bound: 

> Mi(Co)p^, (35) 

for some positive constant pi = pi{Co). 

Case II: Suppose ^ > ^, ^ > ^. By geometric constraint, ^ = 1 + 0{Cq'^). To estimate 
Q3 from below, we decompose the vectors {(^i,^2'~(^i + ^2))} i^^o the orthogonal sum of 
{(^i, — li, 0)} and {(^'1, ^'1, —2^'i)}. Then Q3 is expressed into the sum of three terms as: 



Write: 



^3(^1,^2) = ((^1,^2, -(^1 + ^2)), G3(^i, ^2, -(^1 + ^2))' 
= ((ll,-|l,0),G3(|i,-|i,0f) 

+((^i,i'i,-2^i),G3(^;,$;,-2^;f) 
+ 2((ii,-ii,o),G3(^;,^;,-2^;)^). 




1 1 
1 





V 0/ 



then the bar term of (36): 

^1-10 
((li,0,0); 1 



1 



V(0) - V12 V12 Vi3 
Vi2-V(0) V(0) V23 
V V13-V23 V23 V(0) ) 





V oy 



(36) 



2(li,(V(0)-Vi2)li) >ci/||i|2, 



(37) 
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where c here and after will denote a positive constant depending only on po- Also 1 is a 
shorthand for d x d identity matrix. Similarly, we express: 



and write the prime term by Lemma 2 as: 



v 



1 
1 1 

-2 1 





V 0/ 



/ 



((^i,0,0), 



1 1 

1 




?1 



V(0) + Vi2 - 2Vi3 Vi2 Vi3 

Vi2 + V(0) - 2V23 V(0) V23 

u \ ) \^ Vi3+V23-2V(0) V23 v(o) y V 

= (6V(0) + 2Vi2 - 4Vi3 - 4V23)^i) 

= (I'l, 8(V(0) - Vi3)|'i) + (^'i, (2(Vi2 - V(0)) + 4(Vi3 - V23)) ^'i) 

>c4l^il'-c(/+n24~')l€'ll' 

> c4l€'iP(l - c {{pT-^^f + (prr3^))) > ci4|^i|2. 
The mixed term is equal to : 

' V(0)+Vi2-2Vi3 V12 Vi3 

((Ii,-Ii,0), V12 + V(0) - 2V23 V(0) V23 
^ Vi3 + V23 - 2V(0) V23 V(0) y V y 

= (li, (V(0) + V12 - 2Vi3)^i) - (li, (V12 + V(0) - 2V23)^i), 

= (ll,2(V23-Vi3)0, 



) 



and so is bounded by: 



Thus: 



\mixed term\ < C2\^i\ • • ?'i2?'i3 



C-i 



(38) 



(39) 



Q3 = ^3(^1,^2) > ciP^lliP + ci4|^l|2 - C2|^;| • |li|ri24" 



(40) 



15 



The mixed term may then be controlled by the positive terms through the following Young's 
inequality: 

< ^^-P^l^il^ + ^^^^-^l^'iP, (41) 
with 6 a small number in (0,1). Then, since p/ri^ < 2Cq"^, it follows that for any ^ < 2, 
{p/fizf '^ < for Co large enough. Thus, 

l^'il • \lMs' < \e ■ P^\U + \0 ■ ri,\i',\\ (42) 

which allows the mixed term to be absorbed into the positive bar and prime terms. Combining 
(I^H), we conclude that: 

Qz{^iA2)>¥\li? + crU^'^\ (43) 
which in the original (^1,^2) variables reads: 

^3(^1,^2) > cp^l^i - ^2l' + crUil + ^21'- (44) 

We finish the proof with inequality ( ^4[ ) and ( p5| ) . □ 
Proposition 2, N=4 Case 

We now turn to = 4, for which inequality ( |4^ ) is very helpful. Let r„, n = 1, 2, 3, 4, be four 
distinct points in R"^, d> 2, and assume that ri2 is the minimum length p. Then we show that 
there is a positive constant c depending only on po so that the minimum eigenvalue of G4 is 
bounded from below by cp^. 

Proof: We order and r4 according to the lengths of the three sides intersecting at them. The 
longest length at r4 is larger than that at r^. If they are equal, then the second longest length 
at r4 is larger than its counterpart at r^, and so on. Generically, we are able to order and 
this way. Now rj, i = 1,2,3,4, determine a tetrahedra in R'^. Due to geometric constraint, r23 
and ri3 are on the same order. So are ri4 and r24. With no loss of generality, we can assume 
that ri3 = r23 = a, and ri4 = r24 = /3. Let r34 be 7, which satisfies the inequalities: 

7<a + /3, /3<a + 7;a</3. (45) 
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We consider all the possibilities under (pq). 

Case I. Suppose 2 > ^ > C^^, where C2 > is a large constant to be selected. We have four 
subcases: I 1.1: 1 < ^ < Ci and 1 < < Co; I 1.2: 1 < ^ < Ci and > Cq; I 2.1: ^ > Ci 
and 1 < ^ < C'q- I 2.2: | > Ci and 1 < ^ > C^. Case II: ^ < C^^, which imphes with (H) 
that 1< ^ < We have two subcases: II 1.1: 1< < Co and II 1.2: 1< > Co 

As in the analysis for = 3, we assume that p is smaller than po- The II. 1 is very similar to 
the first case of = 3, in that all lengths are comparable to each other. Writing 



/ 



(^l,l2,€3,-(^l+^2+^3)) 



V 



10 
10 
10 
-1 -1 -1 



(^1, ^2,^3,0)' 



then: 



Q4(^l,^2,^3) = (^l,^2,^3) 
2(V(0) - V14) V(0) + V12 - Vi4 - V24 V(0) + Vi3 - Vi4 - V24 

V(0) + V12 - Vi4 - V24 2(V(0) - V24) V(0) + V23 - V24 - V34 

\ V(0) + Via - Vi4 - V24 V(0) + V23 - V24 - V34 2(V(0) - V34) / 

Using lemma 1 again, we can factor out with remaining matrix being positive and bounded. 
We find that there is n = ^{Cq, Ci, C2) such that: 



\ / . \ 



^2 

V^3; 



Q4 > 



(46) 



Now for I 1.2, we decompose {(^i,^2)^3> + ^2 + ^3))} ™to the orthogonal sum of 



{(Ii,-Ii,0,0)} 



and 



Then: 



{(^'l,^'l,^'2,-2^'l-^2)}- 



Q4(^i, ^2, ^3) = (di, -li, 0, 0), G4(li, -li, 0, 0)^ 
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+ 2((|i, -li, 0, 0), G4(^;, ^1, ^2, -2^1 - ^ 



Writing: 



v 



-^1 







110 
10 
1 








we see that the bar term is equal to: 



(li,2(V(0)-Vi2)li) >cp^|| '2 



ll 5 



Writing: 



^1 

^2 

V -2^i - $2 ; 



1 
1 


-2 



1 



1 1 
-1 1 



^2 



V 0/ 



the mixed term is equal to: 



(ll, -2(Vi4 - V24)^'i) + (ll, (Vi3 - V23 + V24 - Vi4)^^; 



Similarly the prime term is equal to: 
/ 

((^'1,^2), 

\ V23 - 2V24 - 2V34 + 2V(0) 



8V(0) - 8V24 2V23 - 2V24 - 2V34 + 2V(0) 

2V(0) - 2V34 



+ ((^i,^2) 



2Vi2 - 2V(0) + 4(V24 - V14) Vi3 - V23 + V24 - V 



Vi3 - V23 + V24 - Vi4 

The first matrix of (|50|) can be expressed as the product: 
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^ -1 \ / 



1 -1 
1 



V(0) V23 V24 
V23 V(0) V34 
V V24 V34 V(0) J 



I 1 







1 

-1 -1 1 



2 
1 
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hence is positive definite and bounded from below by a positive constant Hi{Ci,C2) times 
^i4l(^'i'^2)P- It follows that: 

^4(^1,^2,^3) > cp^llil' + A^i(Ci,C2)(|^iP + 1^^2)4 - cp4"'llil(ICil + 1^2!) 

-c(/ + pri:'M[f - cp{ri:' + • 1^2!) 

> cp^llil' + C2M[f + > cpi, (51) 

where the mixed term is handled as for iV = 3 with Young's inequality and Cq is chosen large 
enough for given Ci and C2. 

We now consider I 2.1 and I 2.2. Decompose {(^1,^2,^3, ~(^i + ^2 + ^3)} i^^to the orthogonal 
sum of {(|i,|2,-(li +l2),0)} and {(^'1, ^i, ^'1, -3^1)}. Then: 

Q4(^l,^2,^3) = ((ll,l2,-(ll+l2),0),G4(|i,|2,-(ll+l2),0f) 

+(((^i,^'i,^i,-3^'i),G4(€i,€'i,^'i,-3€if) 
+ 2((li,l2,-(li + l2),0),G4(^'i,^'i,^'i,-3^'if). (52) 

Write: 







\ 


/ 


1 











\ 


/ 




\ 













1 












I2 






-(I1+I2) 






-1 


-1 


1 














v 





/ 


V 











1 


/ 


V 





/ 



Then the bar term is equal to: 

.J 2(V(0)-Vi3) V(0)+Vi2-Vi3-V23 \ 

\ V(0) + V12 - Via - V23 2(V(0) - V23) / 

which is larger than: 

c(/|li-|2l' + 4lli+l2l'), (53) 
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by applying and the = 3 result. We express: 



(ri,ri,ri,-3^'i) 



/ \T 



1 











1 


1 








1 





1 





CO 








1 






V oy 



and so: 



V(0) + Vi2 + Vi3 - 3Vi4 Vi2 Vi3 Vi4 

Vi2 + V(0) + V23 - 3V24 V(0) V23 V24 

Vi3 + V23 + V(0) - 3V34 V23 V(0) V34 

\^ Vi4 + V24 + V34 - 3V(0) V24 V34 V(0) J 



The mixed term is equal to: 



^ I -1 ^ 



2(li,^2,0,0) 



1-10 
10 
1 








V(0)+Vi2+Vi3-3Vi4 V12 Vi3 Vi4 

V12 + V(0) + V23 - 3V24 V(0) V23 V24 

Vi3 + V23 + V(0) - 3V34 V23 V(0) V34 

y Vi4 + V24 + V34 - 3V(0) V24 V34 V(0) J 






V oy 



2(ll,l2 



V12 - V23 - 3(Vi4 - V34) V12-V23 I / ^'i 

V12 - Vi3 - 3(V24 - V34) V(0)-V23 )\0 
2(ll, (V12 - V23 - 3(Vi4 - V34))^'i) + 2(^2, (V12 - Vi3 - 3(V24 - V34))^'i), (54) 



which can be written as: 

= 2(li + I2, (V12 - V23 - 3(Vi4 - V34))^'i) + 2(^2, ((V23 - V13) - 3(V24 - Vu))^;). (55) 
It follows that the mixed term is bounded by: 

c|li +I2I • I^'iki3(max(r^2"\r^3"^) +r^4'V(C2)) 
+c|l2ki2?'i3"^|lil +c\^i\ ■ \l2\ri2ri4^- 
20 



The prime term is equal to: 

(^;, (12V(0) + 2Vi2 + 2Vi3 + 2V23 - 6Vi4 - 6V24 - 6V34)^i) 

= (i;, 18(V(0) - Vi4)|'i) + (^'i, (2(Vi2 - V(0)) + 2(Vi3 - V(0)) 

+ 2(V23 - V(0)) + 6(V24 - V34) - 12(V24 - Vi4))^i) 

> c-4i^;p - c-(4 + 4 + 4Mi\' - ^(C2)(r234"' + n24"')l^'il' 
= - /x(C2)cr^ - /x(C2)cri) > c4l^i|^ 

if Ci is chosen large enough for given C2. In case of I 2.1, the mixed terms involving r^J^ 
can be controlled by a Young's inequality as in A?" = 3, using Ci sufficiently large. The 
terms ''i2''i3"^|^'i| • |^2l ^13 '^^{''12"^' ''is'^ll^'il ' l^i +^2! can be estimated by (Co)Pr^2 = 
^12 ■ (.^0)^^12 ' (P ~ 1}), times the ^ bar or prime factors, then using again Young's 

inequality, thanks to the relatively large coefficient r^^ in front of In other words, we use 

Ci being much larger than any chosen Cq. Observe that |^2p ^ |l^2 ~ ^iP + il^2 + ^iPj 
that the mixed terms are again controlled by the prime and bar terms. In case of I 2.2, we make 
Cq itself large to control the term ri2r^^^|^'^| • \^2\- The other terms involving ri4 are standard 
and controlled by large Ci. Note that if C € (0, 1] 

r C-i C-1t c/2 c/2 I ri3^ 2 

ri3max|r^2 ^Hs I = '^is ^^12 — 



Thus when multiplied to + ^2! " l^'il it is bounded by 

2^31^1+^2! +^g—ri2\^l\ < 2^31^1+^21 +2^14l^lh 

with Ci much larger than chosen Cq. If C € (1,2), ri3 max{r^^^, r^^^} = r^^, and its product 
with 111 + I2I • l^'il is bounded by f4||i + I2I' + ^l^'iP < f ^^3!^! + ^2!" + f'^Ll^il', since 
Cf^ < for large Ci. Summarizing the above, we conclude that: 

^4(^1,^2,^3) > + c/lli - 61' + c^lli + ^2l', (56) 
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which, in (^1,^2)^3) variables, is: 



<54(^l, $2, ^3) > C (/l^l - ^2!' + + ^2 - 2^31' + + ^2 + ^3l') 



(57) 



Finally we consider II. The case II 1.1 is no different from I 1.1. Notice that for II 1.2, we 
have essentially two separate scales /3 >> 7, thanks to a and /3 being on the same scale. 
Decompose {(^i,l2>^3> -(^i + ^2 + ^3))} into the orthogonal sum of -|i,l3, -I3)} and 
-^i, -^i)}- The bar term is: 

(di, -li, I3, -I3), G4(li, -li, I3, -Izf)- 



By writing: 



^3 

V } 

we simplify the bar term into: 



1 
-1 1 









1 

-1 1 





I3 



3j> 



2(V(0)-Vi2) V13 + V24-V14-V23 
Vi3 + V24 - Vi4 - V23 2(V(0) - V34) 

Then the bar term is bounded as: 



U3 



(58) 



2(ll, (V(0) - Vi2)|i) + 2(^3, (V(0) - V34)|3) + 2(^1, (Vi3 - V23 + V24 - Vi4)|3) 



> cp^||i|2 + cr3^^||3|2 _ ^IIJ . |v,3 _ Vi4 - (V23 - V24)| 

To obtain the last inequality we used lemma 3: 

l-C/2 l-C/2 

|Vi3 - Vi4 - (V23 - V24)| < cri2r34rL-' = cr^f rg^i^ 



(59) 
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The last term is small for large Ci,C2 when Q < 2. Applying Young's inequality yields the 
same bound as (pg|). 

Next the prime term is simplified by using: 



^1 



10 
110 
-10 10 
-10 1 



?1 








The prime term becomes: 



(^;, (2V(0) + 2Vi2 - 2Vi3 - 2Vi4 - 2V24 - 2V23 + 2V(0) + 2V34) ^'1) 

= (^i, (8V(0) - 2Vi3 - 2Vi4 - 2V23 - 2V24) ^'1) 
- (^;, (4V(0) - 2Vi2 - 2V34) ^'i) > c4ii;|2 



(60) 



The mixed bar-prime term is: 
( \ ( . . . . \ / 







1-10 
0100 
01-1 
0001 



V(0)+Vi2-Vi3-Vi4 V12 Vi3 Vi4 
Vi2+V(0)-V23-V24 V(0) V23 V24 

Vi3 + V23 - V(0) - V34 V23 V(0) V34 

\^ Vi4+V24-V34-V(0) V24 V34 V(0) y 
= (ll, (V23 - Vi3 + V24 - Vm)^;) + (I3, (Vi3 + V23 - Vi4 - V24)^'i) 

Hence the mixed term is bounded by: 

c{n2ri^^ + n2rii^M'i\ ■ llil + c(r34r^J^ + rs^r^^^M'il ■ II3I. 



^1 






(61) 



All the terms in (61) can be estimated as before with Young's inequality, and we have: 



g4(^i,^2,^3) >cp^llil' + c4ie 



3l ' 



which is: 

Q4(li, ^2, ^3) > cp^l^i - ^2l' + 7^12^3 - ^1 - ^2l' + + ^2l'- (62) 

Summarizing all the cases, we finish the proof of the proposition. □ 
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3 Properties of the N-Body Convective Operator 

We now define a sesquilinear form /ijv[*Ar, ^at] for '^n,^n € L'^ ^O®^^, by the expression 

hN[^N,'^N]= [ dK VR^iv(R)-GAr(R)-VR$Ar(R). (63) 

and a quadratic form /ijv[^Ar] = ^Ar[*iV) *iv]- We take as the form domain 

D(/ijv) = {*iv e L'^ (n^^^ : G C°° (n^^^ , supp^jv C 0®^^ for some A;, 

and ^'jv(R)=0 for R G (64) 

Here we made use of an increasing sequence of open subsets of defined as 

n'^^k = {Ke : p{K) > h. (65) 

Clearly, this form can be expressed as /ijv[*iv, ^*iv] = {'^n,'Hn^n) where T^jv is the positive, 
symmetric differential operator 

d 



1 ^ A d 



(66) 



j j'=l ra,m=l 

with Y){J-Ln) = Y){hN)- Our basic object of interest is the self-adjoint (Friedrichs) extension 
Hn of Hn, which corresponds to the operator with Dirichlet b.c. on dQ,^^ . Note that it 
will follow from our discussion below that the same extension Hn also arises if one chooses 
D(7Y7v) = ^"0° (p.'^^^, rather than as above. The main properties of Hat follow from those of 
the form which we now consider. 

The basic properties of the form are contained in: 

Proposition 3 The sesquilinear form /i7v[^iv, *^*iv] enjoys the following. ■ 
(i) Hn is a nonnegative, closable form. 

(a) For a/Z ^'jv G D(/tjv) and for the same constant Cn in Proposition 2, 

/ijv[*iv] >Cn I dK [p(R)]^|VR^'iv(R)|2. (67) 
(iii)For all G D(/ijv) cind for the same constant Cn in Proposition 2, 

hNi'i^N] > Cn ■ i^^^ j^^^ dK [p(R)]-^|*Af (R)l'. (68) 
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Proof of Proposition 3. Ad (i): non-negativity is obvious from the definition Eq.(|63D and 
Proposition That is closable follows from [14|, Theorem VI. 1.27 and its Corollary 



VI. 1.28. Ad (a): This follows directly from the definition Eq.(^) and the variational formula 
for the minimum eigenvalue of GAr(R). Ad (Hi): For the proof of this inequality, we use 
the Lemma 2 of Lewis |10|. That lemma states that, given an open domain A with smooth 



boundary, then for any function g £ H'^{A) such that Ar.(7(R) > for all R G A and for any 
function ip G Cg°(A) (i.e. = on dA), the inequality holds that 

dR I Ar<7(R)||^(R)|2 

< 4 / dR I AR5(R)ri|VR5(R)P|VRV?(R)|2, (69) 

J A 



This is proved by applying Green's first formula and the Cauchy-Schwartz inequality (see [[L0| ). 
Let us take for each integer k > 1 the domain A^ = ^^^k defined as in Eq.(^). If we define 
g{K) = [p(R)]^, then g G i^^(Afc) for each k (and, in fact, g G C^{Ak)). Furthermore, 

Ar 5(R) = 2C(d - iM-R)]'^ > 0, (70) 
for d > 7 (which certainly holds if C > and d >2) and also 

|Vr<7(R)P = 2C2[^(R)]2C-2. (71) 

If £ D(/iAr), then for some k sufficiently large ^ C'g°(Afc), and all the conditions for the 
inequality ( |69[ ) are satisfied. Hence, we find by substitution that 

dK [p(R)]^|VR^Ar(R)|2 
> / „ dK [p(R)]-^|*^(R)|2, (72) 

whenever ^jy S D(/ijv), for (" > and d > 2. If we now use together (ii) and inequality ([72|), 
we obtain (Hi). □ 



Because of item (i) we may now pass to the closed form hjsf (see |14], VI. 1.4). Its properties 
are given in the following Proposition 4: 
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Proposition 4 The sesquilinear form /iAr[^'jv, $Ar] enjoys the following. ■ 
(i) h]\j is a nonnegative, closed form. 

(a) The domain D{hj\f) consists of the Hilhert space obtained by completion of 

Cq° ^r^®^^ in the inner product 

In particular, (n^^'j C D(/iAr). Alternatively, G D(/iAr) iff^^N G L'^ (n^^y its 1st 
distributional derivative satisfies hN[^N] < oo, and 7^ (^^Afjf^siv^^ = for all integer k > 1, 
where 7^. is the trace operator from (n^^ k) into {^dn^^ f] k) ■ 

(Hi) Both the items (ii) and (Hi) of Proposition 3 hold for his[[^ n] and for all ^ n G D(/iAr). 
Furthermore, 

hNV^N] > CnL-^ ■ ^-^^^W^NWh (74) 
also for all ^'tv G D(/iAr). In particular, h]sr is strictly positive. 

Proof of Proposition 4-' (i) '^s immediate. 

(ii) We first prove the statement tliat (q^^^ C D(/iAr). 

To see tliis, we remark ttiat D(/i7v) is dense in Hq (q^^^ for d > 2. In fact, it is well-known 
that in a bounded domain A the set of functions Cg°(A — T), i.e. functions vanishing on F C A 
in addition to A'^, is dense in -ffo(A) if F is a finite union of submanifolds with codimension 
k > 21. This follows from standard density theorems for Sobolev spaces: see Ch.III of Adams 
|§ or Ch.9 of Maz'ja |l6|. The Theorem 3.23 of |l5| states that C^(A - F) is dense in 



i7g(A) iff F is a (2,/)-polar set, when A = R^. However, the same result is true for any open 
domain A. In fact, repeating Adams' argument, if Cq"(A — F) is not dense in IIq{A), then there 
must be a -u G -f^o(^) element T £ IIq\A), the Banach dual, so that T{u) = 1 but 

T(/) = for all / G C^(A - F). However, by |l5|. Theorem 3.10, this T can be identified 
with an element of D'(A) supported on F. Since this can further be canonically identified with 
an element of D'(R'^) supported on F, the set F cannot be (2,/)-polar. The other direction 
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is even simpler. These arguments go back to [|17| . On the other hand, by Theorem 9.2.2 of 
|16| the set F is (2,/)-polar iff its lower //'-capacity vanishes, Cap (r,i/M =0. A convenient 



sufficient condition for zero //'-capacity is that the Hausdorff {D — 2/)-dimensional measure 
of r be finite, 7/^-^/ (p) < ggg Proposition 7.2.3/3 and Theorem 9.4.2 of (This is 

essentially just the converse of the Frostman theorem, due originally to Erdos & Gillis ]l8[.) 
In the case considered, the set F is of Hausdorff dimension D — k, so that ?/^~^' (F) < oo for 
k > 21 {= for k > 21). Thus, the set F has zero //'-capacity as required. Clearly, D(/iAr) 
defined in the statement of the Proposition 3 above coincides with Cq° (^Q^^ — F^ , where the 
set F = {R G f]®^ : r„ = Tm,n ^ m} has codimension = d>2. Therefore, taking D = Nd, 
/ = 1, k = d and A = 17®^ we obtain the density of D(/iAr) in (n^^^ , as claimed. 

As a consequence, for any £ Hq ^0®^^ there exists a sequence of elements S 
D(/ijv) converging in //^-norm to \I'7v- Next, we observe that 

hNl^N] < BnW'^nWhi (75) 

for some coefficient Bj\f > 0. This may be proved by using the variational principle for the 
maximum eigenvalue A™'^^(R) of GAr(R) and then the continuity in R of A™'^(R) over the 



compact set Q"^^ to infer A^^^(R) < B^- This inequality states that the //^-norm is stronger 
than the /lAr-seminorm. Thus, convergence in norm of ^'j^^ G Y){hj\f) to ^'at G Hq (q®^^ 
implies both that ^^pp^ — > in L"^ and also that HnI^^^^ — ^n^] —> as m,n ^ cxd. Com- 
paring with iQ, Section VF1.3 we see that this means precisely that ^'a? G D(/i7v)- Therefore, 
(n'^^^ C D(/iAr). This is the first statement of (ii). 

Next, we recall from Section VL1.3 that D(/i7v) is characterized as the Hilbert space 
obtained by completion of T){h]\j) in the inner-product (|73|). Since this 
is certainly contained in the Hilbert space defined in (ii) above. However, since we have shown 
that (j^®^) C D(/iiv), the completions of D(/iAr) and (o'^^) are the same. 

Finally, we prove the alternative characterization of D(/ijv) in (ii). We note by Proposition 
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3(ii) that for each ^'jv £ D(/ijv) and for each k 

< A;^C^^.||^;v||h^. (76) 

Thus, the i?/i^-norm is stronger than the ^il®^fe^-norm on D(/ijv)|f^0jv^. By definition, for 
each G D(/iAr) there is a sequence ^'^^ G D(/i7v) converging to in -fT/ij^-norm. This 
sequence must also then converge to ^'Tvln^iv in [Q,^^ j.)-norm.. Passing to the hmit in 



(|7^, one then obtains its vahdity for all D(/iAr). This implies that D(/i7v) C ^il®^^ 

for each integer k. Furthermore, the trace 7^ onto the codimension-1 set (9il®^)p|il®^fc 
is continuous from H 

^(n^^k) into H^/^ i^{dn^^)n^^^k)- Since £ D(/i^), we see 

that 7fe ( ) = and, passing to the limit, 7/c (^ArL^iv ) = as an element of 

H^/^ (^{dn®^)r\n^^k)- That is the "only if" part of the characterization. The "if" part is 
very standard. For each ^'jv obeying the alternative set of conditions and A; > 1, we may 
define '^^^^ by extending the restriction ^fTvlQaiv^ again to the whole of O®^, defining it to be 
outside of ^2®^^. Because of the conditions on "^n, the new function ^-jj'^ G (n^^) for 



each k > 1. See Theorems 3.16 and 7.55 of |15|. Thus, ^'^^ G D(/i7v) for all A; > 1. However, 



(1 -Xn«iv^) ^'iv||^ (77) 



where Xcs^j. is the characteristic function of Q'^^ k. Because ||^'Ar||/ijv < 00 by assumption, 
the righthand side goes to zero by dominated convergence as A; — > 00. Thus, we conclude that 
limfc^oo II^AT — ^'^^ll/ij^ = 0, which implies that G D(/iiv)- 



For (Hi): We note that the righthand side of inequalities (|67| ) and ( pS]) in Proposition 3 (ii) 
& (in) are just certain weighted //^-norms and L^-norms, respectively, and both of these are 
bounded by the /ijv-norm on Y){h]\f). Thus, the argument used to extend inequality ( |76|) from 
D(/iAr) to D(/iAr) applies also to extending (|67|)-(|68|). Noting that p(R) < diamO = L for all 
R G n^^, we derive inequality (^ from (|8|). □ 

Remark: The proof does not work for Q, = T^, the d-dimensional torus. In that case, inequality 



(67) of Proposition 3(ii) is still valid, where p(R) = min„_^„ i^gg^ — '''m + L ■ k| has now 
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period L in each direction as required. Unfortunately, the function ^(R) = [p(R)]'' does not 
belong to (^(^d^®N^ away from the set V where /o(R) = 0. It has singularities also on the 
codimension-1 set F' of points where |r„ — r^l = — fml, with r* a periodic image of r„. 
Unless the domain of D(/iAr) is chosen to be = on F', these singularities would contribute a 
surface term in the Green's formula, invalidating (|6^). However, if that condition on D(/i7v) 
is imposed, then the resulting closed form hj^ has Dirichlet b.c. on T', which is unphysical. 
On the other hand, we expect that these are really just problems with the proof and that the 
inequality ( |68[ ) still holds with periodic b.c. Methods used to derive general Hardy-Sobolev 
inequalities ([|^, Ch.2) should apply. 

We now exploit the previous results to study the Friedrichs extension T^tv of Wat. Its 
existence is provided by the First Representation Theorem of forms ( []T^ , Theorem VI. 2.1) 
which states that there is a unique self-adjoint operator TIn whose domain is a core for 

JiN and for which JimI^n.^n] = {'^n,'Hn^n) for every G D(^Af) and <I>Ar G I^{'Hn)- We 
now discuss the essential properties of this operator that we will need later: 

Proposition 5 The Friedrichs extension TIn enjoys the following: 

(i) I-Ln is strictly positive, with lower bound > CnL~'^ ■ ^'^~2^ ■ 

(ii) The spectrum o/TIn is pure point. 



Proof of Proposition 5: Ad (i): (74) and |14], Theorem VI. 2. 6. Ad (ii): We use the Corollary 



to Lemma 1 of Lewis [10|. His hypothesis "HI is satisfied by the increasing sequence ^2®^^ for 
integer fc > 1. His hypothesis 712 is true with = and Ck = Cn ■ as a consequence of 
(|7^. Finally, his third hypothesis holds, with the role of his function p{x) played by Cm-^—t^ ' 



[p(R)] and = C^ ^'^ ■ k by (^). Lewis' proof exploits the Rellich lemma for the 
domain to show that the identity injection / : Hhj^ is compact, by 

approximating it in norm with compact operators /^(^'Ar) = ^'jv^ defined above. The segment 
property holds for $7®^^, since its boundary is C°° except for a finite number of corners where 
the two parts of its boundary, {R E fi®^ : p(R) = k} and (90®^, intersect. □ 
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4 Proofs of the Main Theorems 



We now prove the main results of the paper, using the properties oiTiN proved in the preceding 
section. We start with: 

Proof of Theorem 1: By a stationary weak solution of (|^ at k = 0, we mean a sequence of 
G L2 (o®^) indexed by > 1, such that, for each > 1 and for all e D(?YAr), 

{nN^N,QN) = {^N,Gy) (78) 

where for N >2 

G^(R) = ^7(r„)G^_i(.. .?;...)+ E F{rn,r^)e*N^2{-^-^-) (79) 

™ pairs {nm} 

is the inhomogeneous term of Eq.(^) and G^(ri) = /(ri). Because this quantity for > 1 
involves the correlations of lower order, our construction will proceed inductively. We may 
assume that G (in fact, G*j^ G away from the set T). This statement 

is true for N = 1 and, for N > 2, may be assumed to be true for all M < if the statement 
in Theorem 1 is taken as an induction hypothesis. Only the above regularity property of GJ^ 



will be used in the induction step. Thus, it is enough to show that (78) has a unique solution 



Oat for any Gn G {q.®^^ for each A^ > 1. That is, we must show that for each A^ 

{'Hn^n,Qn) = {^n.Gn) (80) 

for all $Ar G D('HAr) has a unique solution Qn £ L'^ for any chosen Gat. 

We shall first show that Gat G D(/iAr) where &n is any weak solution of ( |80| ) with Gat G L^. 
To do so, we introduce the smoothing operators 

S% = {l + enN)-\ (81) 

In terms of the resolvent operator R{z,A) = {A — z)~^ this may be written as 

S% = -R{--,nN^=R{-l,eHN). (82) 
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These smoothing operators have the following properties: First, they are self-adjoint operators 
with < 1 for all e > 0. Second, because Hn is closed and — ^ is in its resolvent set, 

it follows from the first equality of (|8|) that S% : (vl^^^ ^{'Hn)- This exhibits the 
"smoothing" property of the 5^. Third, 5^ for each e > commutes with Wat, or, more 
correctly, Sj^Ti^ C TCnSJ^. Finally, because lim^^o ^Wtv = in the strong resolvent sense, it 
follows from the second equality of (^2|) that 

lim ||5^^Ar - ^'7v||i2 = (83) 

e— +0 

for all^N & (r?^^). We now observe that, if ©at satisfies ( pO] ) for any Gn in L"^-, then for 
any e D(W7v), 

hN[^N,S%QN] = {Hn^n.S^Qn) 

= {S]^^n,Gn) = {^n,SI,Gn). (84) 

In particular, if we apply this to <I>Ar = Sj^@n, then we find for the quadratic form h]\f[Sf^QN] = 
{Sj^@]\f,<S%G]\f) and, thus, 

hNlSf^ON] < ||©7v||l2 • WGnWl^ (85) 
uniformly in e > 0. Since, in addition, the form h]\f is closed and s — limg^o^^QAf = ®N by 



Eq.(||), it follows from Theorem VI.1.16 of |14] that Qn e D(/iAr). 



In that case, for any G Y){Hn), the equation (|^) may be rewritten 

hN[<^N,&N] = {<^N,GN). (86) 

Furthermore, T>{TCn) is a core for D(/ijv) by the First Representation Theorem for forms: see 
|14|, Theorem VI. 2.1, item (ii). By the same Theorem VI. 2.1, item (in), it follows from ( p6| ) 



that ©TV G D('H7v) and that 

HnQn = Gn (87) 
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with equality as elements of L (f]®^) . We observe, since 'H^ is bounded, that the equation 
( |87D is equivalent to 

Gat = Hj^Gn (88) 

However, for any Gat G {ti®^^ the righthand side of (^) exists, again by boundedness of 
and it defines an element Qn = T~(-]^^Gn G ^0~(-n)- Thus, the weak solution exists and is 
unique. □ 

Proof of Theorem 2 (i): The proof of existence and uniqueness here very closely parallels the 
previous one, but is even easier. For this reason, we will discuss only a few details. As in the 
previous case, we may begin by introducing a symmetric sesquilinear form, 

N 

in6 



N 

h^;;^\^N,^N]=hN[^N,<^N] + Y. dR (-ArJP/2^'^(R) • (-Ar„r/2cI>^(R). (89) 



densely defined on either the same domain as before, D (^h^j^^^^ = D(/i7v), or, with identical 
results, D (/iS^"^) = C§° (r?®^). Clearly, this is the same as /ij^''^ [^jv, ^iv] = (^jv, wj^'^^^iv), 
where ^iv"'' is the differential operator in Eq.(D with D (h^^"^^ = D (/ij^''-'). We may now 

~(k ) 

consider the self-adjoint (Friedrichs) extensions of these operators, denoted "H^y^ , just as before. 
We may observe that there is a basic inequality, 

with some constant Aj\[ > 0, for all ^'tv S D (^h^^^''^^ . This plays the same role in the present 
proof as inequality ( [7^ of Proposition 4 (in) did in the previous one. It is proved first for 
^(^p)|^^^j g Y) l^/i^^^^, by expanding the elements of D (h^^^'^^ in a series of eigen- 

functions of the Dirichlet Laplacian (— Ar)^), which are complete in Hq(^Q®^^. Then, the 
result is extended to h^^''^"^^] by taking limits. Note that the inequality (|90|), in particular, 
implies that the operator TL^^''^ is strictly positive, with lower bound > KpA]\f / L'^p . Therefore, 



the inverse operator '^jy*'^ is bounded, as before, and unique weak solutions Q^^'"" of the 
stationary equations are easily constructed with its aid. 



-1 
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A last point which requires some explanation is the regularity 0^''''* G Hq (q^^^ of solu- 

(k )* )\ 

tions. In fact, it follows as before that Qjy'' £ D ihjy'' ). It therefore suffices to show that 

We may identify D (h^^^^^ as the completion of the pre-Hilbert space 
with the inner product 



(^JV,^7v),(.p) 



N 



(91) 



See |14], Section VI. 1.3. However, we have the elementary inequality 



- N 


p/2 


- N 


.n=l 


< Cn,p 


.n=l 



(92) 



with Cn,p = N^'P ^^/^ for p > 2 and = 1 for I < p < 2. Using then the Parseval's equality 
for Fourier integrals, it follows that the norm H^'a' || , (Kp) is stronger on Cq° (Q®^] than the 
Sobolev norm 



Since Hq (q^^^ is defined to be the completion of Cq° (q'^^^ in the norm 
that D (h'^N"^) C H[; , as required. □ 

Proof of Theorem 2 (ii): To construct the weak-L^ limits of 0^''''* for Kp 



(93) 

\hv, it follows 
► 0, the main 



thing that is required are a priori estimates on the L^-norms uniform in k„ > 0. These are 



provided as follows. First, we note that 0^''''* G D(/iAr) because Q'^^^'" G 
D (n^^"^) C HP (p.^^^ C D(/iiv) for p>l. Thus, we may apply Proposition 4 (iii), inequality 



and 



(74), to calculate that 



0(«;p)*||2 
, N IIl2 



IN 



— '^N^ ^'^N I 



(94) 



with C'j^ = [CN{d - 7)2/2]"^ In other words, 

\\e[";'''^*\\r2 < c'mL'^\\g[7'^*\ 



(95) 
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Using the expression (79) for G^'' in terms of the lower-order Qj^f , for M < N, it fohows 
that 

l|giv^^^*llL^ < N • 11711^.(0) ||ei;^\* 11^. + \\Fh^inmMN'^2\\L^- (96) 

It is then straightforward to prove inductively from (|95|) and (|9|) the main L^-estimates 

\\G^N''^*\\LHn^^) < 2 • i2KNFLY'' ■ N\ (97) 

and 

\\QN'^*\\LHm^) < {2KnFLY ■ ^! (98) 

where Kn = maxi<Af<iv C^,/ and F = max{||/||i2(Q), ||-^||/2(s^^n)}- 

A further crucial estimate may be extracted from the preceding discussion. Using the 



inequality /iat B^''^* < ||6)v'''''||l2 ||G)(^*''''||j;^2 contained in (p^) and the bound on G)^ , 
(PW), it follows that 



,(Kp)*|| ||^(Kp)*ll • /|n7j|\ T .1 „ t2 u^„^a r'i'^p)* 



hN 



< 



{2KNFL^f-W] . (99) 



In other words, we have a uniform bound for the 0^*"^* in the norm of the Hilbert space H^j^: 

\\e^;^'''*\\h^ <2-{2KNFL^f -m (100) 
This follows by combining estimates (|9l) and (p9|). 



We now consider the weak-L^ limits of 0^''''* as Kp — > 0. We note first, because of the a 
priori bound ( p8| ) and weak compactness of the unit ball in L^, that any sequence Kp""* — > 
contains a subsequence Kp" ■* such that w — lim^/^oo ©jv'' ^* ~ ©tv^* exists, and 

W&TWL^n^n < {2KnFLY ■ N\ (101) 

Furthermore, because of the additional a priori estimate ( |100| ), we may extract a further sub- 
sequence Kp which converges weakly in Hh^^, and the limit then satisfies 

\\e^S^*\\h^<2-{2KNFL^f -Nl (102) 
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We wish to characterize ah the possible such weak sequential limits Q)y . If the weak hmits 

(n') 

along subsequences hip diiTQ all identical, then, in fact, the weak limit exists and equals the 
unique subsequential limit. 

We shall show that, in fact, all of the weak subsequential limits coincide with 0-^, the 
unique weak solution of the zero-diffusivity problem. First of all, we observe that for all > 1 
and all G D (t^S^^^ 



N 

;;^-)*), (103) 



because Q^^^^* is a weak solution of the pt/i-hyperdiffusivity equation. As a consequence of The- 
orem 2(i), we may take D (h^n"^^ = D (/iS^^^) = (jl®^). Because D(wS^^^) C D(wS^''^), 
( |103| ) is therefore true for all G (Q'^^V independent of the value of > 0. Passing to 



the limit along subsequence Kp" , we then obtain 



{■HN^N,eT) = {^N,GP*), (104) 

for all <I>Ar S ^il®^^. This is not quite the statement that 0^'** is a weak solution of the 
zero-diffusivity equation, with our definitions. For that to be true it is required that ( 1041) hold 



for all <I>7v G D('HAr). By the same argument as above, Cq° (^Q,'^^^ is a dense subset of D(/iAr) 



in the Hilbert space Hhj^- Thus, we would like to take the limit in Hhj^ to obtain (104) for all 



$Ar S Y)(Ti.N), as required. To do so, however, requires that ej^^* G D(/i7v), so that we may 
write 

^(0)* 



{'^N,G'Z), (105) 



In this form, the limit may be taken to obtain (|104D for all <I>Ar S Y){TIn)- Thus, to complete 
the proof, it is enough to show that 0^^* S D(/iAr). 

(0)* 

To demonstrate the latter regularity of ©jy , we shall use the second characterization of 
D(/iAr) in Proposition 4:(ii). We already have the estimate ( [1021 ). All that is required in addition 
is to show that 

35 



for all A; > 1. To obtain this, we remark that for each k the identity injection : H^^ 
H'^ (i}^^ is compact for any s < 1, because the identity injection from H^j^ (^^^^^^ to 
(n^^k) is bounded by (^) and the identity injection (l^®^fc) into H' (n^^ k) is com- 
pact, by the Rellich lemma. We may use the above compact embedding for each fixed k to 

(n'") 

extract by a diagonal argument a further subsequence Kp such that 



lim 

n'"— >oo 



= (107) 



for all k > I. However, for each k, the trace 7^ is continuous as a map from into 
L2 (aO®^n^^®^fc) when s > 1/2. Furthermore, 



7. 1 e'/'"'" 



(108) 



for all n'" . Thus, passing to the limit, we obtain ( |106| ). □ 

5 Concluding Remarks 

We make here just a few remarks on some further results of our analysis and some outstanding 

problems for future work. 

(i) Regularity of the Solutions 

The construction above produces solutions 0^ G L away from the 

singular set T. In fact, as was mentioned in the Introduction, it is expected that 0^ are Holder 

regular, Q*j^ S (q^^^ . Such additional regularity of the solutions of the singular-elliptic 

equations may follow from Harnack inequalities 20 1. 

(a) N -Dependence of Spectral Gap and Invariant Measure on Scalar Fields 

The Proposition 2 has only been fully proved here for < 4. Assuming that it holds for general 

N, the question of the A^-dependence of the constant Cn appearing in its statement has also 

some importance. As we have seen, the solutions 0^ constructed for k = obey an L^-bound 

WQliWmn^^) < B"" ■ N\ (109) 
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in which B is proportional to the inverse of minjv>i Cat. If Cn is bounded from below uniformly 
in N, then the above constant B < oo. In that case, the correlation functions determine a 
characteristic functional via the series 




(110) 



absolutely convergent for ||'!/'IIl2(q) < -B ^. A measure fi* on scalar fields 9 G L^(J7) such that 

<^*{i^) = I e'^^'^^fi*{d9), (111) 

JL2(n) 

is therefore uniquely determined by the correlation functions. That such a measure actually 
exists is a consequence of the Minlos-Sazonov theorem (see [^], Theorem V.5.1), if it can 
be shown, for example, that <I>* defined by Eq.( [L10D is a positive-definite, weakly continuous 
functional on L'^{Q) and ©2 is the kernel of a trace-class operator on L'^(Q), i.e. Jq dr ©2(r, r) < 
oo. The latter would follow from the regularity discussed in (i). 

The measure fj,* so constructed would be the natural candidate for an invariant measure 
on the scalar fields. Whether the dynamical equation Eq.(Q) itself can make sense for k = 
with Dirichlet b.c, even in a suitable weak sense, is an unresolved issue. However, for a 
periodic domain, or = T*^, the d-dimensional torus, there should be a sensible theory of weak 
solutions to Eq.(|^) and we conjecture that the reconstructed measure ^* will be invariant under 
realizations evolving according to that equation. 
(Hi) Time- Dependent Solutions and Relaxation to the Steady-State 

From our work in this paper there follow some further results for time-dependent solutions of 
the correlation equations, Eq.(^). In fact, by standard semigroup theory (see ||T^, Ch.IX), a 
unique solution to Eq.(^) with initial datum 0Ar(O) G ^$7*^^^ may be (inductively) obtained 
via the Riemann integrals 

ONit) = e-*^^eiv(0) + /* ds e-(*-^)^^G7v(s) (112) 

Jo 

with G]\f given by Eq.([79[), in terms of the strongly continuous, contraction semigroups Tjv(t) = 
g-tHjv_ -yy-g i-gfpain from precise theorem statements here. Furthermore, because of the strict 
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positivity of spectrum oi TCn established here, the semigroup T]\f{t) is strictly contractive and 
the limit exists 

^lim \\@Nit) - ©|vr||L2(Q<g.jv) = 0. (113) 

Thus, the time-dependent solutions converge strongly in to the stationary solutions con- 
structed in this work. All of the results on existence of solutions for Kp > and their conver- 
gence to zero-diffusivity solutions for Hp 0, which were proved above for stationary solutions, 
also carry over to the time-dependent solutions. 
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